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We report on a search for the process B + -> K*°(892)K + using 232 x 10 6 r(4S) -> BB decays 
collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. From a 
signal yield of 25 ± 13 [stat] ± 7 [syst] 5+ -> K*°(892)(-» K~tv+)K + events, we place an upper 
limit on the branching fraction B(B + -> K*°(892)K + ) of 1.1 x 10~ 6 , at the 90% confidence level. 
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PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



We present a measurement of the branching fraction 
B(B+ -> K* a {892)K+) based exclusively on B+ decays 
to the final state K + K~ir + . Charge conjugate states 
are assumed throughout. In the Standard Model (SM), 
B — > K* K decays are dominated by b — > dss gluonic 
penguin diagrams (see Figure 1 (a) in [1] ; for the charged 
decay the spectator d is replaced with u). Such tran- 
sitions provide a valuable tool with which to test the 
quark- flavor sector of the SM (see, for example, [2-4]). 
The mode B+ -> K*°(892)K+ is also relevant for the 
interpretation of the time dependent CP asymmetry ob- 
tained with the B° — > (j)K® mode. To leading order the 
CP asymmetry equals sin2/3 for this mode, where (3 is the 
Unitarity Triangle (UT) angle. However, sub-dominant 
amplitudes proportional to V* b V us could produce a devi- 
ation AS^, K o from sin2/3. Exploiting SU(3) flavor sym- 
metry and combining measured rates for relevant b — > s 
and b -> d processes (including B+ -> (892)i^+) , a 
method is introduced in [5] to place a bound on AS^ro . 
Measurements yielding a significant deviation in excess 
of such a bound would be a strong indication of physics 
beyond the SM. Furthermore, B+ -> K*°(892)K+ is one 
of several charmless decays that can be used, together 
with U-spin symmetry, to extract the UT angle 7 [6]. 

Theoretical predictions for B(B+ -» K*°(892)K+) in- 
clude B(B+ -> K*°(892)K+) > (0.46±8;g|) x lO" 6 [2] 
and B(B+ -» i?* (892) J ftT+) w 0.5 x 10~ 6 [3]— both us- 
ing SU(3) flavor symmetry and experimental information 
for charmless B decays, and -» K*°(892)K+) w 

0.31 x 10~ 6 using perturbative QCD factorization [4]. 
Prior to the analysis presented here, the only experimen- 
tal limit placed on B(B+ -> )?* (892).fC + ) was that pre- 
sented by the CLEO collaboration at the 90% confidence 
level (CL) [7]: B{B+ -» i?*°(892)if +) < 5.3 x lO" 6 . 

The data used in this analysis were collected with 
the BABAR detector [8] at the PEP-II asymmetric-energy 
e + e~ storage ring at the Stanford Linear Accelerator 
Center. Charged particle trajectories are measured by 
a five-layer double-sided silicon vertex tracker and a 40- 
layer drift chamber located within a 1.5 T axial magnetic 
field. Charged hadrons arc identified by combining en- 
ergy loss information from tracking (dE/dx) with the 
measurements from a ring-imaging Cherenkov detector. 
Photons and electrons are detected by a CsI(Tl) crystal 
electromagnetic calorimeter. The magnet's flux return is 
instrumented for muon and neutral hadron identification. 

The data sample consists of (232 ± 3) x 10 6 BB pairs 
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collected at the T(45) resonance (on-resonance data), 
corresponding to an integrated luminosity of 211 fb -1 . 
It is assumed that neutral and charged B meson pairs 
are produced in equal numbers [9]. In addition, 22fb _1 
of data collected 40 MeV below the T(45) resonance (off- 
resonance data) are used for background studies. 

B meson candidates are reconstructed from three 
charged tracks. The charged tracks are required to have 
at least 12 hits in the drift chamber and a transverse 
momentum greater than 0.1 GcV/c. They arc fitted to a 
common vertex; momentum must be conserved at this 
vertex. Two of the tracks must have opposite charge 
and a signal in the tracking and Cherenkov detectors 
that is consistent with that of a kaon. We remove tracks 
that pass electron selection criteria based on dE/dx and 
calorimeter information. 

We perform full detector Monte Carlo (MC) simula- 
tions equivalent to 2.4 x 10 5 signal B+ -> F*°(892)(^ 
K~ir + )K + decays. For background studies l.Oab -1 of 
generic BB decays are simulated, as are over 100 exclu- 
sive B meson decay modes (~ 10 4 — 10 6 events/mode), 
approximately half of which are charmless. MC samples 
are generated with EvtGen [10], while the detector re- 
sponse is simulated with GEANT4 [11]. All simulated 
events are reconstructed in the same manner as data. 
Off-resonance data are used to measure the properties of 
the light quark continuum decays, e + e~ — > qq (q = u, d, 
s, c). 

For correctly reconstructed signal events, 
AE = E* B — \fs/2 peaks at zero, while toes = 
\J (s/2 + po • pb) 2 /Eq — p 2 B peaks at the B mass. The 
resolutions of these largely uncorrelated kinematic vari- 
ables, for signal events, arc sa 20 MeV and ss 2.5 MeV/ c , 
respectively. Eg is the B meson candidate energy in 
the center-of-mass (CM) frame, Eq and y/s are the 
total energies of the e + e _ system in the laboratory 
and CM frames, respectively, and po and ps are the 
three-momenta of the e + e~ system and the B meson 
candidate in the laboratory frame. The distributions 
for continuum events are slowly varying. For B meson 
decays in which particle misidentification occurs, the AE 
peak is shifted by 500 or more MeV. Events are selected 
with 5.22 < toes < 5.29GeV/c 2 and \AE\ < 0.1 GeV 
The AE restriction helps to remove background from 
two- and four-body B meson decays at a small cost to 
signal efficiency. 

Continuum quark-antiquark production is the domi- 
nant background. To suppress it, we select only those 
events where the angle 9t in the CM frame between the 
thrust axis of the B meson candidate and the thrust axis 
of the rest of the event satisfies |cos6>t| < 0.9. For 
continuum events, which tend to be jet-like in the CM 
frame, the distribution of | cos 6t \ is strongly peaked to- 
ward unity whereas the distribution is uniform for sig- 
nal events in which little kinetic energy is available in 
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the CM frame. The number of continuum background 
events present per signal event is reduced by a factor of 
approximately two with the application of the | cos 9t \ 
cut. We also construct a Fisher discriminant [12] J 7 , a 
linear combination of five variables: the zeroth and sec- 
ond angular moments of the energy flow — excluding the 
B candidate — about the B thrust axis; the absolute value 
of the cosine of the angle between the momentum vector 
of the reconstructed B candidate and the beam direc- 
tion; the absolute value of the cosine of the angle between 
the thrust axis of the reconstructed B candidate and the 
beam direction; and the output of a multivariate, non- 
linear B meson candidate flavor tagging algorithm [13]. 
The Fisher coefficients are obtained from samples of off- 
resonance data and B + — > K + K~ir + MC. A loose cut 
of \T\ < 3.0 is applied. This cut eliminates a negligible 
fraction of signal and background events and is applied 
only to define a range for the fit (see below). 

Further discrimination between signal and continuum 
background is achieved by utilizing the variables mj^ 
and cos8h- The invariant mass of the K*°(892) candi- 
date, niK-w, is restricted to 0.744 < my, < 1.048 GeV/c 2 . 
We also require that the cosine of the helicity angle, 9h , is 
less than 0.9, where 9h is defined to be the angle between 
the pion track and the spectator kaon track in the rest 
frame of the S*°(892). This angle depends on the spin 
of the intermediate resonance: for B + — > K + K ~ir + via 
the spin-1 S*°(892) resonance, the distribution of its co- 
sine is quadratic. For continuum events the distribution 
is approximately uniform. At high cos Oh the final state 
pion has low momentum and is difficult to reconstruct. 
This causes a sharp drop-off in the efficiency between 0.9 
and 1.0. The selection criterion cos6># < 0.9 makes an 
unbinned fit to the variable possible at a cost of losing 
5% of signal events. 

After the selection described above, the B + — > 
A 7 *°(892)H K-tt+)K+ selection efficiency is 26%. In 
signal MC studies, the signal candidate is correctly re- 
constructed 94% of the time. The remaining candidates 
come from self-cross-feed (SCF) events that stem from 
swapping one or more tracks from the true B meson de- 
cay with tracks from the rest of the event. 

To identify backgrounds from B meson decays the se- 
lection criteria described above are applied to the MC 
samples. Using the efficiencies of the selection criteria 
and world average branching fractions we find that the 
largest expected contributions derive from b — > c transi- 
tions and from charmless 3-body decays in which a kaon 
is misidentified as a pion or vice versa. The b — > c events, 
combinatoric in nature, are continuum-like in most of 
the fit variables. As such, the contribution to the fitted 
signal is small: <~ 1% of the number of b — > c events 
expected to be present. For charmless 3-body sources, 
however, the contribution to the signal yield — as a pro- 
portion of the number of events present — is considerably 
larger. This is particularly true of B + — > K + K~K + and 
B+ -> pK+. The decay B+ -> K+R-R+ includes non- 
resonant and several intermediate resonance states, in- 



cluding the narrow (f> state. In order to increase statistical 
precision, <j>K + is considered separately from the rest of 
the K + K~K + final state, which is modeled using the re- 
sults of the Dalitz plot analysis in [14]. Since the <f> is very 
narrow, its interference with the other K + K~K + states 
can be neglected in this context. The modes B + —> <pK + , 
B+ -> K+K-K+, and B+ -> pK+ , arc therefore in- 
cluded as components of the fit. This eliminates biases 
on the fitted signal yield due to these channels. 

The contribution to the signal yield due to all other 
sources of B meson background (including the b — > c 
modes discussed above and numerous charmless modes) 
is estimated from simulation at 3 events — to which a 
conservative ±100% uncertainty is assigned. The uncer- 
tainty accounts for poorly known branching fractions and 
simulation limitations. Individually these sources con- 
tribute at a low level. As such, rather than including 
components for each of them in the fit, a correction to 
the fitted signal yield is made. 

It is also necessary to consider backgrounds from B 
meson decays that have the same final state as the sig- 
nal mode. MC studies of the B+ -> K + K-n+ Dalitz 
plot show that the only contribution that needs to be ac- 
counted for is B+ -> K* °(U30)(-> K--k+)K+. Wc use 
the LASS parameterization for the Sq°(1430) lineshapc, 
which consists of the X5°(1430) resonance together with 
an effective range non-resonant component [15]. We take 
B(K^{U30) -» K-n + ) to be equal to § x (93±10)% [16]. 
A maximum likelihood fit to three variables — toes; AS, 
and T — is performed in a region of the K~ir + invari- 
ant mass spectrum between 1.048 and 1.800 GeV/c 2 in 
an analogous way to how we fit the main signal (see be- 
low). A 90%_CL upper limit of 2.2 x 10~ 6 is placed 
on B(B+ -> S*, (1430)S+)._From simulation and the 
branching fraction B(B+ — > K 5°(1430)i^ + ) as obtained 
above, we estimate — assuming zero interference — that 
5tl 4 B+ -» if* (1430) H K-tt+)K+ events will be 
present in the region 0.744 < vtik-k < 1.048 GeV/c 2 , con- 
tributing 2±2 events to the fitted B+ -> S*°(892)(-> 
K~n + )K + signal yield. The central value of this esti- 
mated contribution is calculated from the central value 
of B(B+ -> Kf(U30)K+) while the uncertainty on the 
contribution covers the contribution obtained from the 
upper limit on B(B + — > Sq°(1430)S + ) and uncertain- 
ties in the parameterization of the S5°(1430) lineshapc. 

A correction is applied to the fitted B+ -» _S*°(892)(^ 
K~-K + )K+ yield to account for B+ -> if*, (1430) H 
K~ti + )K + ; including a component in the fit is ineffec- 
tual since B+ -> S^°(1430)(-> K-n+)K+ is signal-like 
in the majority of the fit variables. 

An unbinned extended maximum likelihood fit to the 
five variables rn.ES, AS, T, vcik-k-, an d cos##, is used 
to extract the total number of B+ -> S*°(892)(^ 
K~ir + )K + and continuum background events. The like- 
lihood for the selected sample is given by the product of 
the probability density functions (PDFs) for each indi- 
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TABLE I: PDFs are described by one or more of the follow- 
ing functions: ARGUS [17] (A), Breit-Wigner (BW), Crystal 
Ball [18] (CB), double Gaussian (DG), exp(P) + C— where P 
is a polynomial in the fit variable and C is a real scalar (E), 
Gaussian (G), linear (L), one-dimensional non-parametric [19] 
(N), quadratic (Q), Voigtian — a Gaussian convolved with a 
Breit-Wigner (V). 



Component 


m E s 


PDF variable 
AE T m Kl < 


cos 9h 


Signal 












Truth-matched 


CB 


DG 


DG 


BW 


E 


SCF 


CB 


L 


DG 


G+L 


E 


Background 












Continuum 


A 


L 


DG 


BW+L 


Q 


B+ -> 4>K+ 


CB 


G+L 


DG 


V+E 


E 


B+ -> K+R-K+ 


A+G 


G+L 


DG 


Q 


N 


B+ -> pK+ 


CB 


G+L 


G 


G+L 


E 



vidual candidate, multiplied by the Poisson factor: 
1 , N 

1=1 

where N and N' are the number of observed and expected 
events, respectively. The PDF Vi for a given event i is 
the sum of the signal (S) and background (B) terms: 

V t = ((1 - /) (V s ) i + / {V§ OF )J + J>f (*f )< > 

3 = 1 

(2) 

where N s and are the yields for the signal compo- 
nent and the background components j, and / = 0.06 is 
the fraction of SCF signal events (treating true signal and 
SCF separately reduces the correlation between toes an d 
cos Oh for signal from 14% to 1%). The four background 
terms comprise the continuum distribution and the three 
B meson background modes described above. The PDF 
for each component is the product of the PDFs for each of 
the fit input variables: V = V mES Vi\ E V r V mK ^ V CO se H - 
Any correlations between the variables are such that bi- 
ases brought about in the fit are negligible and we treat 
each PDF as independent and uncorrelated. 

The PDF forms are presented in Table I. The parame- 
ters of the signal and B meson background PDFs are held 
fixed to the MC values. The parameters of the contin- 
uum PDFs are allowed to float except for the endpoint of 
the ARGUS function. The signal and continuum yields 
are floated in the fit while the three B meson background 
yields are fixed to their MC expectations. 

Individual contributions to the systematic uncertainty 
are summarized in Table II. The systematic uncertainties 
that arise from fixing PDF parameters for the signal and 
B meson background components are estimated by vary- 
ing these parameters, one at a time. Correlated param- 
eters in the relevant PDF are adjusted accordingly and 



the maximum likelihood fit is repeated with the shift in 
the signal yield taken to be the systematic uncertainty. 
The parameters are varied cither by the lcr uncertain- 
ties obtained when evaluating them from MC or such 
that we account for any discrepancy observed between 
data and MC (whichever is larger). Such discrepancies 
are identified with the calibration channel B + — > D°ir + , 
which has a topology similar to the signal and a much 
higher branching fraction. For V s os g H , rather than vary- 
ing the PDF parameters, a second-order polynomial is 
used. This is the expected shape when neglecting detec- 
tor and reconstruction effects. For the non-parametric 
PDF, its smoothness is varied. The positive and negative 
shifts for each varied PDF parameter/shape are added 
separately in quadrature. The same procedure is used for 
the fixed yields of the B meson background modes and 
for the SCF fraction, /, which is varied by ±20% (rel- 
ative). The uncertainty due to fixing the ARGUS end- 
point for the continuum background toes PDF is found 
by floating this parameter and observing the shift in the 
signal yield; this shift is found to be less than a hun- 
dredth of an event. The systematic uncertainties associ- 
ated with the subtraction of events from the fitted signal 
yield due to B+ -> Kf(U30)(--> K~tv + )K+ and other 
non-K + K ~7r+-final-state B meson decays have been dis- 
cussed above. 

The remaining systematic uncertainty on the signal 
yield is due to possible interference between final states. 
Several thousand MC datasets are produced each con- 
taining B+ i?*°(892)(^ K-tt+)K+ and B+ 
ifo°(1430)(— > K~ir + )K + events that are generated ac- 
cording to their Breit-Wigner and LASS lineshapes, re- 
spectively. Interference between the two modes at the 
amplitude level is modeled. The relative magnitudes 
and phases of the two contributing amplitudes are varied 
randomly between datasets, but the numbers of events 
present in the regions 0.744 < niK-n < 1.048 GeV/c 2 
and 1.048 < m Kll < 1.800 GeV/c 2 are the same for 
each dataset and are equal to the numbers we observe 
in data. The fractional systematic uncertainty is taken 
to be twice the standard deviation of the distribution 
of the fraction fsg2 divided by its average and is found 
to be 8%. For a generated dataset, / 8 g 2 is the modu- 
lus squared of the amplitude of the B+ -> K*° (892) (-» 
K ~ir + )K + mode integrated over the full Dalitz plot, di- 
vided by the modulus squared of the sum of the am- 
plitudes^ the B+ -> j?* (892)(^ K-tt+)K+ and 
B+ -> K* Q (U30)(-> K~-k+)K + modes integrated over 
the full Dalitz plot. 

The uncertainties on reconstruction and selection cri- 
teria efficiencies are evaluated by comparing efficiencies 
for MC and data control samples. A systematic uncer- 
tainty of ±1.4% per track added linearly is taken for 
the particle identification efficiency. A systematic un- 
certainty of ±0.8% on the tracking efficiency is applied 
for each charged track, added linearly. The systematic 
uncertainty on the efficiency of the selection criteria is 
found to be ±13%. The systematic uncertainty on the 
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TABLE II: Breakdown of systematic uncertainties. 



Systematic effect 



Uncertainty 



+ 1.8 



events 



-1.1 
±0.7 events 

+0.4 



-0.6 



events 



Yield 

Fixed PDF parameters 
Fixed SCF fraction 

Fixed B meson background yields in fit 
B meson background contribution 

Non-if+iT"7r + -fmal-state 

K + K _ 7r + -fmal-state 
Final state interference 



±3.2 events 



events 



+2.0 

-5.8 



±2.0 events 



Reconstruction and selection criteria efficiency (RSC) 
Tracking ±2.4% 
Particle identification ±4.2% 
Selection criteria ±f3.0% 



Yield total 
RSC total 

Total number of B events 



-4.7 
-7.1 



events 



±13.8% 
±1.1% 



Total systematic uncertainty (xlO 6 ) on 
B(B+ -> K*°{892)K + ) 



±0.2 



40 
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H 
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total number of B events is ±1.1%. 

A total of 38, 690 events are fitted. The numbers of 
B+ -> <j>K+, B+ -> K+K-R+, and B+ -> pK+ events 
expected in this sample, estimated from simulation, are 
21 ± 2, 23 ± 2, and 4 ± 1, respectively. The yields for 
these B meson background components are fixed at the 
central values. The raw signal yield extracted from the 
fit is 30 ± 13 [stat] ± 3 [syst] events, of which 5^ [syst] 
are estimated to be B meson background events. The 
number of true signal events present in the on-resonance 
data sample is therefore 25 ± 13 [stat] ± 7 [syst]. The 
statistical significance of the result in the absence of sys- 
tematic uncertainties, defined as the square root of the 
difference between the value of — 21n£ for zero signal 
events and at its minimum, is 3.1a. Accounting for sys- 
tematic uncertainties, this significance is reduced to 1.6a. 
The number of signal events is divided by the product of 
the signal efficiency and the total number of B events 
to give the branching fraction B(B+ -> K*°(892)K+) = 
(0.6±0.3 [stat] ±0.2 [syst]) x 10~ 6 . Since the signal is not 
significant, we place an upper limit on this measurement 
at the 90% CL: B(B+ -> K*°(892)K+) < 1.1 x 10~ 6 . 
The likelihood function defined in Eq. (1) is modified 
to incorporate systematic uncertainties through convo- 
lution with a bifurcated Gaussian whose standard devi- 
ations are set to the (asymmetric) total systematic un- 
certainties described above. The 90% CL upper limit is 
then defined to be the value of the branching fraction 
B(B + — > K*°(892)K + ) (which corresponds to a particu- 
lar value of N s ) below which lies 90% of the total integral 
of the modified likelihood function in the positive branch- 



FIG. 1: (a)-(e): distributions of tties, AE, T, tuktv, and 
cos8h. The points with uncertainties show the data. The 
curves show projections of the maximum likelihood fit. A se- 
lection requirement on the likelihood ratio has been applied 
as described in the text. The black, solid curve (no filling) 
shows the sum of all fitted components. The curve with gray 
filling shows the sum of all background components. The 
curve with black filling shows the signal component; (f): like- 
lihood (modified to account for systematic uncertainties, and 
normalized) as a function of B(B + -» K*°(892)K + ). The 
shaded area represents 90% of the total area under the curve 
in the positive branching fraction region. 



ing fraction region. This is illustrated in Fig. 1 (bottom 
row, right). 

The results of the fit to B+ -> K*°(892)(-^ K-n+)K+ 
are illustrated in Fig. 1. The plots are enhanced in signal 
by selecting only those events that exceed an optimized 
threshold for the likelihood ratio R= N s V s /(N s V s ± 
YT j=1 NfVf) where N are the central values of the 
yields and V are the PDFs with the projected variable 
integrated out. 

In conclusion, we have reduced the 90% CL upper 
limit on the branching fraction for the decay B + — > 
iT*°(892)i^+ from 5.3 x 10~ 6 [7] to 1.1 x 10~ 6 . The cen- 
tral value has been measured to be (0.6 ± 0.3 [stat] ± 
0.2 [syst]) x 10 -6 with a significance of 1.6a, consistent 
with the predictions of [2] , [3] , and [4] . This measurement 
can be used to determine an upper bound on AS^o. 
The technique described in Sec. VI of [1] is used, together 
with Eq. (29) of [5] and the findings of the B+ -> 07r+ 
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and B + — > <^>i^T + analyses described in [20] and [14, 21- 
23], respectively, to place a 90% CL upper bound of 0.11 
on |AS , 0^o|. Systematic uncertainties on the branching 
fractions used to determine this bound are accounted for. 
The bound presented here is significantly more restrictive 
than the bounds of w 0.4 found in [1, 24]. We note that 
these latter bounds are based on fewer theoretical as- 
sumptions (see [5]). We have also placed an upper limit 
of 2.2 x 10~ 6 on the previously unmeasured branching 
fraction of the decay B+ -> K* °(U30)K+. 
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